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Abstract 

We study spontaneous symmetry breaking in ((p'^)i^i using the hght-front for- 
mulation of the field theory. Since the physical vacuum is always the same as the 
perturbative vacuum in light-front field theory the fields must develop a vacuum 
expectation value through the zero- mode components of the field. We solve the 
nonlinear operator equation for the zero-mode in the one-mode approximation. 
We find that spontaneous symmetry breaking occurs at Acriticai = 47r (3 + Vs) , 
which is consistent with the value Acriticai = 54.27 obtained in the equal time 
theory. We calculate the value of the vacuum expectation value as a function 
of the coupling constant in the broken phase both numerically and analytically 
using the 6 expansion. We find two equivalent broken phases. Finally we show 
that the energy levels of the system have the expected behavior within the 
broken phase. 



I. Introduction 

It was Dirac [l|] who first recognized that different generators of the Poincare 
group could be used as a Hamiltonian for the purpose of quantizing a field theory. 
He showed that, among these, light-front quantization was unique. In light-front 
quantization the Hamiltonian is p~ ' = (/ - p^) 7^2, and p+ = + p^) / ^2 is the 
third or "longitudinal" component of momentum. The longitudinal momentum has a 
positive semi-definite spectrum, and massive excitations cannot mix with the vacuum. 
Therefore, the bare Fock space vacuum is an eigenstate of the full Hamiltonian. 
Since theories such as QCD exhibit spontaneous symmetry breaking, spontaneous 
symmetry breaking must appear through some other mechanism in the context of 
light-front field theory. We will see that spontaneous symmetry breaking occurs in 
light-front theory because the field includes a zero-mode that is not an independent 
degree of freedom. This mode is a complicated operator-valued function of all other 
modes in the theory and may lead to a non-zero vacuum expectation value (VEV). 

In order to investigate the zero-mode in light-front field theory, we will consider 
a discretized 0^ field theory in two space-time dimensions. The formal constraint 
equation obtained from the Dirac-Bergmann quantization procedure relates the zero- 
modes to all the other modes in the problem [Qj. This equation is, however, most 
easily obtained by integrating the equation of motion 0. It is clearly very difficult 
to solve and, to date, only approximate qualitative solutions have been generated 0]. 
The procedure that generates the constraint equation does not specify an operator 
ordering. We will use a symmetric operator-ordering prescription 0. In order to 
render the problem tractable, we will truncate the Fock space to include only the 
lowest-energy mode. We will see that for weak coupling the theory has only the 
trivial solution for the zero-mode. However, as we increase the coupling we reach a 
critical coupling where a pair of nontrivial solutions to the constraint equation appear. 

Section II presents a simple derivation of the zero-mode constraint equation in 
the classical case. Section HI briefiy discusses quantization and mass renormalization 
of the theory. In Section IV we discuss the asymptotic behavior of the zero-mode 
in the large-particle-number sectors. From the asymptotic behavior of the constraint 
equation we show that the theory has a critical coupling. In Section V we use the 
5 expansion to study the solution branches away from the critical point. Section VI 
present some numerical solutions to the equation. We show that the theory has a 
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critical point predicted by the asymptotic behavior of the constraint equation. The 
two nontrivial solutions are plotted as functions of the coupling and the VEV along 
with the 6 expansion solution. We also study the behavior of matrix elements of 
the zero-mode near the critical curves and the energy eigenvalue as a function of the 
coupling. Finally in Section VI we discuss our results and the remaining work that 
is needed on this problem. 

II. The Classical Case 

The details of the Dirac-Bergmann prescription and its application to the system 
considered in this paper are discussed elsewhere in the literature 0, In terms 
of light- front coordinates = (a;° ±x^) /\/2. For a classical field the La- 
grangian is 

£ = - - . (2.1) 

We put the system in a box of length d and impose periodic boundary conditions. 
For most of our discussion we work in momentum space. We define by 

m = A^j2^n{x^y'-''~, (2.2) 

Vd n 

where = 27m /d and the summations run over all integers unless otherwise noted. 
Next, we introduce some notation and separate out the zero-mode. We define 

= — E ^iife • • • Qin ^ii+i2+-+in,0 ■ (2-3) 

11,12, ■■■,inT^O 

Using the Dirac-Bergmann prescription, one can find the canonical Hamiltonian 

Following the Dirac-Bergmann prescription, we identify first-class constraints which 
define the conjugate momenta 

= Pn - ik:^q-n , (2.5) 
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and a secondary constraint, 



n 2 I ^^^O I '^90^2 AE3 



(2.6) 



which determines the "zero-mode" go- This result can also be obtained by integrating 
the equations of motion. One can calculate the Dirac brackets between the coordinates 



d 



d^qkqiSk+i,n 
A k,i 



m 

and with the conjugate momenta, 



[Qm ) Pv 



, m,n ^ 



(2.7) 
(2.8) 



(2.9) 



The total longitudinal momentum is given by 



(2.10) 



One can show that P'^ has vanishing Dirac brackets with qq, and P : 



. 



(2.11) 



III. Canonical Quantization 

To quantize the system one replaces the classical fields with corresponding field 
operators. One uses commutators instead of Dirac brackets and inserts a factor of i. 
One must choose a regularization and an operator-ordering prescription in order to 
make the system well-defined. 

We begin by defining creation and annihilation operators al and Ofe, 

= \ -TTn ' = «-fe > ky^O , (3.1) 
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which satisfy the usual commutation relations 



[ttk, ai] = , 



al, a] 



0, 



(3.2) 



We also define 



An 



(3.3) 



Very general arguments suggest that the Hamiltonian should be symmetric ordered 
[0. This choice of operator ordering produces tadpoles which we eliminate by adding 
an overall constant and a mass counterterm to the Hamiltonian: 



Xd 



1287r2 



^1 

i>o ^ 



Xd 
64^ 



^ 1 ^ O'kO'-k 



l>0 



(3.4) 



Then, the quantum Hamiltonian is 



^ O'kO'lO'm.O'nSk+l+m+nfl 



k,l,m,n 

2 



imn 



Note that the constraint equation for the zero-mode is obtained by taking a derivative 
of P with respect to oq. Consequently, it is natural in the quantum case to order 
symmetric order the constraint equation: 

= ^foo + + 2aoS2 + 2S2ao + 2^ r h 6S3 , (3.6) 

fc>0 ^ 

where g = 2AiTfi^/X. This equation implies that in matrix form is real and sym- 
metric. Moreover, it is block diagonal in states with equal eigenvalues. Using the 
constraint equation, we define a rescaled Hamiltonian H: 
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H 



9677" 



p- 



gj:2 + 9^4 - — — 



+ 



1 ^ ajaoakCii + ajakaotti 



+k+l,0 



12 2 2 

1 ^ aka^a^k + d-kO'o^k — o-o 
4^ k 



(3.7) 



IV. One Mode, Many Particles 

It is reasonable to assume that the lowest-energy mode will be the most important 
one in the constraint equation (3.6). We therefore study the case where the zero-mode 
is just a function of the lowest-energy mode. In this case, the zero-mode is diagonal 
and can be written as 

ao = fo\0){0\ + J2fk\k){k\ . (4.1) 

fc>0 

Equivalently, one can think of the zero-mode as an operator-valued function of the 
number operator = a^a. The VEV is given by 

(O|0|O) = -i=(O|ao|O) = -i=/o. (4.2) 



Substituting (4.1) into the constraint equation (3.6) and sandwiching the con- 
straint equation between Fock states, we get a recursion relation for /„: 

= 9fn + fn + (4n - l)fn + {u + I) fn+1 + ufn-l . (4.3) 

If we take /o = and we assume that n/„_i evaluated at n = vanishes P], then 
we see that all the /^'s are zero. This corresponds to the unbroken phase, = 0. 
Our objective is to determine whether spontaneous symmetry breaking occurs and a 
nonzero solution for oq appears as we increase A (decrease g). 

We begin our analysis of Eq. (4.3) by finding its asymptotic behavior for large n. 
If /n ^ 1 in this limit, then the term will dominate and 
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from which we deduce that 



(4.4) 



J.im /„ ~ (-l)"exp(3"constant) . (4.5) 

We now argue that we reject this rapidly growing solution. The part of the 
Hamiltonian (3.7) corresponding to the lowest mode is diagonal in and the zero- 
modes will affect all of the energy levels. If /„ is large for large n then the high 
energy levels will be strongly affected. The paradigm for spontaneous symmetry 
breaking is the symmetric double well which has a ground state localized in either 
well rather than at the symmetry point. This paradigm indicates that the behavior of 
the system is unaffected by the barrier for energies far above the barrier separating the 
wells. Hence, we only seek solutions where /„ is small for large n. This is the central 
condition that will be used to determine the critical couplings in all the subsequent 
calculations. We therefore neglect the term for large n; it is the terms linear in /„ 
that dominate, giving 

/n+l+4/„ + /n-l = 0. (4.6) 

There are two solutions to this equation: 

/„oc(v^±2)". (4.7) 



We reject the plus solution because it grows with n. Dropping the cubic term from 
(4.6) we define the generating function 

oo 

F{z) = Y: fuz"" . (4.8) 

n=0 

If /„ goes like then the radius of convergence of F{z) is 2 + ^3 and we 

expect F{z) to be singular at \z\ — 2 + y/S. Similarly, if /„ ~ (-\/3 + 2^ , then we 
expect F{z) to be singular at \z\ — 2 — -\/3. 
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The function F{z) satisfies the differential equation 

1 dF{z) _ g-l + z 
F{z) dz ~ ~ z^ + Az+1 

whose solution is 



m ~ V 2-^3 J I 2 + V3 j 



(4.9) 



(4.10) 



Note that this solution for F{z) has singularities at the expected values of z. If 
we want /„ to have the asymptotic behavior — 2^ for large n, then we must 
eliminate the branch point of F{z) at \z\ = 2 — y/S. This gives the condition 

'^~^^^=K, K = 0,l,2. (4.11) 



2^3 

Only K = gives g > 0. Therefore, we find a critical coupling 

S'critical = 3 - V3 . (4.12) 

We conclude that there is a critical value of the coupling constant for which there 
is a nonzero value for ao and a solution to the linearized equation for /„ exists that 
does not grow rapidly for large n. We can compare our critical value of g with that 
obtained in the equal-time formulation. Chang ^ finds that Acriticai = 54.3 which 
differs from our result 



Acriticai = 47r (3 + v^) ^ 59.5 (4. 13) 

by about 10%. 

Of course, we need to determine if there is just an isolated critical point or if there 
is a continuous range of values of g < (^critical for which Oq has nontrivial solutions. 
This requires that we investigate the full nonlinear equation. Away from = where 
the ttg term can make substantial contributions we will use both the 5-expansion and 
numerical methods to answer this question. We denote the values of /o vs. g that 
satisfy (4.3) the "critical curve" . 
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V. The ^-expansion 

The 5-expansion is a powerful perturbative technique for hnearizing nonhnear 
problems. It has been shown to be an accurate technique for solving problems in 
differential equations, quantum mechanics, and quantum field theory [j^. 

We rewrite Eq. (4.3) as 

- 1 + 4n) U + fn^'' + {n + l) + nU^^ = 0. (5.1) 

Setting 6 = gives the linear finite difference equation which is the zeroth-order 
approximation in the 5-expansion. One then expands in powers of 6 about 6 = 0. 
One recovers the problem of interest aX 6 = 1. Expanding about 5 = we have 

g = g(^) + 6g('^ + ... , /„ = /(o)+5/(i) + ... , n = l,2,... (5.2) 

and 

f'n^'' = /„,(l + 51n/„^ + ...) 

= /f +5(/r + /rin/rV... . (5.3) 

Substituting this into equation (4.1) we find, to zeroth order in 6, 

(^?(°) + 4n) /f + in + l) + n/f 1 = . (5.4) 

This zeroth-order equation is the same equation as that studied in Section IV with 
g displaced by one. As discussed before, we can determine the solution for large n. 
One then finds that 

/r = /o(v^-2)" , g^'^ = 2-V3. (5.5) 

To first order in 6 we obtain an inhomogeneous second-order finite difference 
equation: 

{g^'^ + 4n) /J^) + {n + l) + = (in + . (5.6) 
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This equation can also be solved exactly. We obtain 



1 In f 2 - v^l 



X < 



2" (2 + V3) 

(2 + V3) E^--^ 

p=l 



2p 



n 

E- 

p=iP 



/o 

2V^' 



(5.7) 



The second term in (5.7) grows with n so we demand that its coefficient vanishes, 
which gives 

^(1) = -ln/2 + -Lin (2 + VS) (2 - VS) , (5.8) 



/J^) = -;^(x/3-2)"ln(2 + x/3) 



(5.9) 



These results are plotted in Fig. 1. The dashed line is a plot of the critical curve for 
5 = 1, 



g^(2-Vs) (|l + -Lln(2 + v^)j -ln/o% 



(5.10) 



away from /q — 0. The expansion behaves badly near fo — because of the Iu/q 
term in the expansion. The (5-expansion analysis clearly shows that there is a critical 
curve and not merely a critical point. 



VI. Numerical Solution 

We can study this critical curve in detail by looking for numerical solutions to 
Eq. (4.3). The method used here is to write Eq. (4.3) as a set of M simultaneous 
equations: 



= (^7-l)/o + /o' + /i, 

= (^ + 3)/i + /f + 2/2 + /o, 

= (^ + 7)/2 + /| + 3/3 + 2/i, 

= (5 + 4M-1)/m + /m' + M/m-i . 
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(6.1) 



In the Mth equation /m+i is set to zero. Since we seek a solution where is 
decreasing with n, this is a good approximation. We then pick a value of g and look 
for real solutions for /o, /i, . . . , /m- We find that ioi g > 3 — VS the only real solution 
is /n = for all n. For g less than 3 — \/3 there are two additional solutions and near 
the critical point |/o| is small and 

^ /o (2 - 73)" (6.2) 

As g decreases (A increases), the solution for |/o| increases. The critical curve is 
indicated by the sohd line in Fig. 1. The solution of (6.1) converges quite rapidly 
with M. The critical curve is approximately parabolic in shape: 

^ = 3 - - 0.9177/o^ (6.3) 

For a given value of /o and g Eq. (4.2) can be used to calculate all values of /„. 

It is interesting to study the behavior of the constraint equation (4.3) away from 
the critical curve. In Fig. 2 we plot |/„| as a function of n and /o for g = 1.2. We 
see that, as n becomes large, all the increase and as /o approaches the critical 
curve, which is at /o ~ 0.2700 for g = 1.2, all the l/nl's decrease rapidly. As /o 
increases beyond the critical curve the |/„|'s increase rapidly once again. The fact 
that \fn\ increases rapidly on both sides of the critical curve is a manifestation of the 
nonlinearity in (4.3). 

We can also study the eigenvalues of the Hamiltonian (3.7) for this one- mode 
problem. The Hamiltonian is diagonal in the number operator so the energy 
eigenstates are just the eigenstates of N. Thus, 

(n| H In) = ln{n - 1) + - | - + + ^fl, . (6.4) 

In Fig. 3, the dashed lines show the first few eigenvalues as a function of g without 
the zero-mode. Observe that the vacuum is at zero for all g. When we include the 
zero-mode, the energy levels shift as shown by the solid curves. The vacuum is at zero 
energy for g > (^critical but at g = (^critical there is a phase transition and the energy 
decreases below zero as g is decreased. We also see that for g < ^fcriticai all the higher 
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energy level increase above the value they had without the zero-mode. The higher 
levels change very little, as our paradigm would suggest, because /„ is small for large n. 

VII. Discussion 

In the context of (t)\j^i field theory on the light front, our paradigm for spon- 
taneous symmetry breaking suggests that spontaneous symmetry breaking occurs 
when fields can develop a zero-mode. This zero-mode gives rise to a nonzero VEV 
and the full vacuum remains the perturbative vacuum. In the broken phase, the 
theory behaves exactly as expected and the numerical value of the critical coupling 
Acriticai = 47r (3 + v^) agrees well with the value obtained in the equal-time theory 



Spontaneous symmetry breaking is simpler to understand in light-front field theory 
because the entire effect comes from one mode. However, the problem of solving for 
this one mode is quite difficult. In the literature, it has been suggested that a direct 
solution of the zero-mode problem may be intractable f^. We hope to have convinced 
the reader that this is not true. 

We interpret the existence of more than one solution to the constraint equation to 
be spontaneous symmetry breaking. However, this spontaneous symmetry breaking 
is non- dynamical. We have no motivation to choose one solution to the constraint 
equation over the others. The conventional argument that when the system is coupled 
to a heat bath it will dynamically pick the lowest energy state to be the vacuum does 
not apply here. 

Many problems remain to be addressed. A more complete solution of the zero- 
mode problem including all the oscillators would be very interesting, if only to confirm 
our results. In reference the authors retain the one particle and, implicitly, some 
of the two particle states for all of the modes and find a solution for the critical cou- 
pling Acriticai = 47r (3.184. . .). Also, some questions regarding operator ordering still 
remain. When we chose a quantum Hamiltonian, we demanded that it be symmetri- 
cally ordered and we treated the zero-mode as an ordinary field operator. However, 
go is not an ordinary field operator and can, in principle, be written in terms of the 
other field operators oq = co + Y.CkO^kCik + ■ ■ ■■ In this sense, the Hamiltonian we chose 
is not really symmetrically ordered after all. It is unclear whether this is a problem. 



'Critical 



54.27 i, |T0 . 
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Figure Captions 

Figure 1. g — 247r//^/A vs. /o = ViTrVEV. The solid curve is the critical curve 
obtained from numerical solution of (6.1) with M — 10. The dashed curve is 
the critical curve obtained from the first-order 5-expansion. 

Figure 2. as a function of n and /o for g = 1.2 from the numerical solution of 
(6.1) with M = 10. 

Figure 3. The lowest three energy eigenvalues as a function of g from the numerical 
solution of (6.1) with M = 10. The dashed hue is the symmetric solution fo — 
and the solid line is the solution with /o 7^ for < gicriticai- 
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